Understanding the saturation of proton-driven Weibel instabilities and implications 

for astrophysics 
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The linear growth rate and saturation level of magnetic fields for Weibel instabilities driven by 
ion temperature anisotropy, defined as a — (Tx/Tm) — 1 where T± and Tji are ion temperatures 
perpendicular and parallel to the wave vector, are derived in the small a-limit. ft is shown that the 
ratio of the saturated magnetic energy to the initial ion energy scales as the fourth power of the 
electron to ion mass ratio, m/M, for an initially unmagnetized plasma with a < M/m. Particle- 
in-cell simulations confirm the mass scaling and also show that the electron energy gain is of the 
same order of magnitude as the magnetic field energy. This implies that the Weibel instabilities 
cannot provide a faster-than-Coulomb collisionless mechanism to equilibrate ion-electron plasmas 
with ions initially much hotter than electrons, a key component in low-luminosity astrophysical 
accretion flows. The results here also show that the large a-limit formulae used in the study of 
magnetic field generation in collisionless shocks are only valid if a > M/m. 



I. INTRODUCTION 
A. Weibel instabilities and astrophysical shocks 

Weibel instabilities, whereby an initially anisotropic 
plasma distribution relaxes to convert the free energy in 
the anisotropy into magnetic fields e.g. [H, 0, H| have re- 
ceived much recent attention in the study of astrophys- 
ical plasmas, particularly in the context of influence on 
the observed emission associated with collisionless astro- 
physical shocks. These include shocks of the crab pulsar 
wind [4j, gamma-ray bursts (GRB)3, and between or 
within galaxy clusters 0, [||. The former two are rela- 
tivistic. The roles that the Weibel instabilities may play 
via their ability to grow magnetic fields are (i) to magnet- 
ically trap particles near the shock, facilitating the shock 
formation itself; (ii) to account for synchrotron or jitter 
radiation that would otherwise be absent and hard to ex- 
plain without a source of field growth; (iii) to suppress 
thermal conductivity, as needed to explain observed tem- 
perature inhomogeneities Q in the ion-electron plasmas 
of Galaxy clusters; and (iv) possibly producing a faster- 
than Coulomb coupling or electron acceleration in the 
post-shock regions near shocks within or between Galaxy 
clusters. This is important for the estimate of Galaxy 
cluster masses and therefore fractional dark matter mass 
fraction estimates The latter are inferred from mea- 
sured electron temperatures, the assumption that pro- 
ton temperatures are equal to the X-ray inferred electron 
temperatures, and virial equilibrium. 

The plasma compositions for the astrophysical appli- 
cations range from ion-electron to positron- electron pair 
plasmas, to a mixture between the two. In shocks, the 
systematic velocity anisotropy is typically assumed to be 



strong, and is produced by fast particle streaming in the 
shock propagation direction compared to the directon 
parallel to the shock. Current fluctuations parallel to the 
direction of shock propagation lead to current filaments 
that amplify the magnetic field in the shock plane, which 
then further enhances the current filaments. 



B. Weibel instabilities and two-temperature 
plasmas 

We introduce the question of whether the Weibel in- 
stability may also facilitate a faster-than-Coulomb cou- 
pling between ions and electrons in two-temperature plas- 
mas even unshocked. The free energy in anisotropic need 
not be the result of particle streaming perpendicular to 
shocks, but would also result from anisotropy in the local 
distribution of random particle velocities of the species 
supplying the free energy in the absence of shocks. The 
study of the Weibel instability without shocks also has 
a long history [2j but has not been studied for plasmas 
containing ions and electrons at very different average 
eneriges,the motivation for which we now describe fur- 
ther. 

Left alone, a two-temperature ion-electron plasma, 
with the ions at a higher temperature than the electrons, 
will eventually approach an equilibrium in which the 
two species have the same temperature. Coulomb colli- 
sions provide a minimum equilibration rate, but whether 
a faster-than-Coulomb collisionless equilibration mecha- 
nism exists @, H, Q for specific particle distributions and 
external conditions is both a fundamental question in 
plasma physics and important for understading aspects 
of collisionless astrophysical phenomena such as radia- 
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tively inefficient accretion flows (RIAF) (e.g. Ref . [10] ) . 

The survival of a two-temperature plasma for a 
Coloumb equilibration time scale is important for geo- 
metrically thick RIAFs of which advection dominated ac- 
cretion flows (ADAFs) have been the archetype [n], Ell 
[H, [T3]- In contrast to standard geometrically thin ac- 
cretion discs (e.g. Ref. [HI), which can comfortably allow 
10% of the binding energy of the accreting material to be 
converted into observable radiation around black holes, 
two-temperature accretion disks can in principle reduce 
the photon luminosity by many orders of magnitude for 
the same accretion rate. Indeed, the measured luminosi- 
ties from the central engines of some nearby elliptical 
galaxies seem to be 3 to 5 orders of magnitude smaller 
than expected for their estimated accretion rates. Similar 
quiescient phases of accretion are also observed in some 
states of X-ray binaries. There is also X-ray evidence 
that in some supernova remnants the ion temperature is 
much higher than the electron temperature [l6(. 

In the simplest ADAF model, the reduced luminosity 
results from three key features: (1) energy dissipated by 
viscosity is assumed to heat the protons, (2) only elec- 
trons radiate efficiently, and (3) if the only means of en- 
ergy transfer between electrons and protons is Coulomb 
collisions, then a sufficiently collisonless plasma can be 
accreted onto the central object before the electrons re- 
ceive enough energy to radiate the energy dissipated by 
the accretion. The gravitational binding energy of the 
accreting material is then retained as internal energy of 
the hot protons, and can be quiescently advected through 
the event horizon of a black hole. Since the associated 
internal energy per ion in an ADAF is of order m p c 2 ~ 
lGeV, the electrons are required to stay at a tempera- 
ture much lower than the GeV ion temperature to ensure 
a low luminosity for a given accretion rate. ADAFs are 
thick discs because the thermal energy of the ions puffs 
up the disk. More general models of two temperature 
accretion involving outflows and convective feedback are 
also PQDular[17l ll8L Il9j. 

The assumption that heat transfer from the ions to 
the electrons occurs only by Coloumb collisions is neces- 
sary to ensure that a negligible fraction of the dissipated 
thermal energy is transferred to the electrons during the 
accretion infall time. This assumption has yet to be vali- 
dated or invalidated (e.g. Ref. Doing so requires un- 
derstanding the subtle plasma physics of the interactions 
among the ions, electrons, and electromagnetic fields. A 
systematic way to make progress is to test this assump- 
tion under a range of specific conditions. If the assump- 
tion is found to be violated for particular circumstances, 
then the associated instability and faster-than-Coulomb 
coupling would limit the applicability of two-temperature 
accretion paradigm. If the assumption is not violated for 
a particular set of conditions, this would not prove the 
viablity of two-temperature plasmas in all circumstances 
and more cirumstances would need to be tested. While 
the latter issue has discouraged some people from work- 
ing on this problem, we think that a series of rigorous 



tests using different initial and boundary conditions is 
needed: Were the assumption to emerge as consistently 
valid for a range of different initial conditions, this would 
gradually solidify the robustness of the assumption for an 
increasingly wider range of viable conditions, and justify 
two-temperature accretion flows for those regimes tested. 

In this spirit, we study one test of the Coulomb equi- 
libration assumption here, namely, whether ion-electron 
energy transfer can be sped up via Weibel instabilities 
driven by an ion distribution anisotropy. For a real disk, 
the microphysics by which gravitational potential energy 
is converted to ion thermal energy is not fully under- 
stood, but presumably involves turbulence generated by 
a collisionless version 2l| of a magneto-rotatioanl insta- 
bility [2^ . MHD turbulence is believed to anisotropi- 
cally cascade [23| and it is likely in a collisionless system, 
the proton temperature may also incur anisotropies on 
small scales, transiently induced by the cascading tur- 
bulence. Similar to an anisotropic electron distribution, 
an anisotropic ion distribution can also drive a Weibel in- 
stability [HIH which amplifies fluctuating magnetic fields. 
How fast and how much energy the electromagnetic fields 
and the electrons then gain from the isotropizing protons 
is the subject of our present work. 



C. Growth and Saturation of the Weibel Instability 

Since electrons gain energy from electromagnetic fields 
generated the Weibel instability, the electron energy gain 
is closely related to the energy level of the magnetic field, 
the dominant electromagnetic field at saturation. The 
strength of the saturated magnetic field level from the 
Weibel instabilities is therefore important for all of the 
astrophysical applications discussed above: two temper- 
ature accretion plasmas as well as the origin of magnetic 
fields in Gamma-ray bursts [!, H3] and galaxy clusters 
The saturation mechanism is also a fundamental plasma 
physics question. 

If energy transfered to electrons were a significant frac- 
tion of the initial ion energy during an instability growth 
time for the two temperature plasmas discussed in the 
previous section, this transfer would be a much more 
efficient process than Coulomb collisions for tempera- 
ture equilibration. The is because the ratio of maximum 
growth rate of the ion-driven Weibel instability j max (see 
EqJT] below) to collisional ion-electron equilibration rate 
Vei is 
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Here a > is a measure of ion temperature anisotropy 
(defined more precisely in Sec|H|, n the electron density, 
Ad the Debye length, In A the Coulomb logarithm, T e and 
Tj the electron and ion temperatures (multiplied by the 
Boltzmann constant, as throughout the paper), respec- 
tively, m and M the electron and ion masses, respec- 
tively, and c the velocity of light. For a proton-electron 
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plasma with n = 10 12 cm -3 ,Tj = 10 12 if, a = 1, and 
T e = Tjm/M, "fmax/vei ~ 3 x 10 7 . This ratio is large, 
but the question is what fraction of the ion energy will 
be transfered to the electrons when the instability satu- 
rates? 

There have been two saturation mechanisms discussed 
in the literature: (i) the deviation from the equilibrium 
orbits of the particles with an anisotropic distribution 
@; EL HI1 and (ii) the magnetization of the ions Q or 
electrons [(| . These two mechanisms gave the same mag- 
netic field level in the electron-driven case (for either 
ion-electron or pair plasmas) but differ in the ion-driven 
case [13] ■ In the present paper we will show and reaffirm 
that the orbit deviation is the correct magnetic satura- 
tion mechanism for the Weibel instabilities. 

We will derive the growth rate and magnetic saturation 
level in the small anisotropy limit in Sec[Tl] which shows 
that the fraction of the ion energy transfered to the mag- 
netic field via the Weibel instabilities scales as (m/M) 4 
for a single mode. Connection will also be made to previ- 
ous results in the large anisotropy limit. Particlc-in-Cell 
(PIC) simulations in Sec lIIII with reduced ion mass will 
confirm this mass scaling and also will show that the 
electron energy gain is of the same order of magnitude 
as the magnetic field energy. For an initially unmagne- 
tized plasma, this implies that the Weibel instabilities 
are insignificant to the ion-electron temperature equi- 
libration process. However, in Sec lIVI we discuss how 
the transfered energy fraction via the Weibel instabili- 
ties may greatly increase for an initially more strongly 
magnetized plasma. We conclude in Sec. IIV1 



II. CALUCLATING GROWTH AND 
SATURATION OF PROTON-DRIVEN WEIBEL 
INSTABILITIES 

A. The growth rate 

The linear growth rate 7 of the Weibel instability has 
been derived for various anisotropic electron or ion dis- 
tributions 0, IH. For the ion driven case, Ref. [|J derived 
7 in the large anisotropy (growth rate) and low plasma 
temperature limit for the "ion-pinch" case, where the 
electron distribution is an isotropic Maxwellian and the 
ion distribution consists of two bi-Maxwellian counter- 
streams ■ Here we also study instability driven by the 
ions, but we do not invoke any ion bulk streaming, and 
instead allow the instability to develop from ion temper- 
ature anisotropy alone. Also in contrast to (2j, we focus 
on the the small anisotropy (growth rate) , high ion tem- 
perature limit. 

We assume that the electrons have an initially isotropic 
electron distribution given by 



while the ions take on a bi-Maxwellian distribution, 
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Here, n is the density of the electrons or the ions (singly 
charged ions are assumed for simplicity), T e the electron 



temperature, v A 



and T Ix and 2>j_ = T Iy = T Iz 



are the ion temperatures in the x-direction and in the di- 
rections perpendicular to x respectively. To study the 
linear stability of this system against the Weibel insta- 
bility, we assume that the mode wave vector is in the x- 
direction, k = kx, and all perturbations have spatial and 
temporal dependence of exp(ikx — iuot). For simplicity, 
we do not consider any coupling between transverse and 
longitudinal modes [23, [26| and assume that the modes 
are transverse, E = E y y and B = £>z|- 

We follow the procedure of Ref. [2j in linearizing the 
Vlasov-Maxwell system of equations, which leads to the 
following dispersion relation, 

^ + f[i + iz'K.)] + ^[i + i Z '(f I )] + g^'« I )-i = 

(4) 

Here u 2 e = 4wq 2 n/m and oj 2 j = 4irq 2 n/M are the plasma 



frequency of each species, £ e 



u)/(ky/2T e /m) and 



£r = ujj (ky/2Ti x /M), Z' is the derivative of plasma dis- 
persion function Z{£/) = tt^ 1 / 2 exp{— A 2 )/(A — £,)d\ 
[27l ] . and a = Tj±/Tj x — 1 represents the ion tem- 
perature anisotropy. Eqf4] can also be directly derived 
from Eq.(37) in Ref. by removing bulk streaming from 
both species and assuming an isotropic electron distri- 
bution. Verified a posteriori, the growth rate for the 
ion-driven Weibel instability for small to moderate a is 
small, |w|/fcc ~ m/M. Therefore, we take the small ar- 
gument limit, £ << 1, for Z'(£),Z'{£) m -2 - 2i^£, 
to further simplify the dispersion relation Eqf4j In this 
small growth rate limit, the displacement current term 
(the constant 1 on the left-hand side) can also be ne- 
glected. This leads to 
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where (3 = y (1 + a ) is usually much less than 1 

when T e < Ti x and ay/ m/M << 1. Eq. shows that 
the mode is unstable when a > and 



k 2 <k 2 = uLa/c 



(G) 



In fact, this expression for the range of unstable fc's is 
exact and can be obtained from EqU without using the 
small £ approximation of Z . The maximum growth rate 
is reached at k* = /c /\/3, and is given by 
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Eq. (with (3 w 0) was used to derive Eq. (I}. 

From Eq. ((7|) it is clear that the "f m ax is smaller than 
tOpi by a factor of m/M, when T e /m ~ Tj x /M, in this 
small £/ limit, valid for small to moderate a. (For T e ~ 
Tj x , this factor is y/m/M.) Since the linear growth rate 
determines the saturation level of the magnetic field (see 
Sec lII B[) . it is interesting to see how large 7 can become 
in the large £/ limit, valid for large a. Since Z'(£) « l/£ 2 
when £ >> 1, we obtain from Eq|4"{"" 



'Jpi 



T, 



VI + a(kc/u} p i) 



Mc 2 



y/1 + (fcc/wp/) 2 + M/m 



(8) 



which is valid for fe << fco- Using EqJS]we can find that it 
takes a significantly large anisotropy, a >> \/2M/m, to 
achieve £/ >> 1. For large fe, i.e., kc/cjpj » yj M/m, 7 
approaches the maximum value of u> p i y/T I± /Mc 2 , which 
is still less than uj p ] for any non-relativistic Tj±. A rela- 
tivistic calculation using a water-bag distribution found 
that 7 ~ u> r r [24j . Our discussion here indicates that 
this is only possible for relativistic temperatures and ex- 
tremely high a's. For a ~ M/m, EqJT] scales the same as 
the maximum growth rate from Eq[8l 
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It is thus reasonable to infer that Eq[7| has the correct 
scaling for a up to M/m. 



B. Saturation 

The PIC simulations of Ref.[2| showed that the lin- 
ear growth for a particular fc-mode saturates when the 
magnetic field amplifies to B s defined in the following 
expression 



lk ~ y/W±qB s /Mc. 



(10) 



Here, jk is the linear growth rate for the mode k and V± 
and M are the characteristic perpendicular (to k) velocity 
and mass of the instability-driving species (the ions in our 
case) . Reference 0] identifies the right-hand side of EafTUl 
as the magnetic bounce frequency and concluded that 
the saturation results from particle-trapping. Detailed 
analytical and numerical analysis in Ref.[4| showed that 
EqllOl is equivalent to the condition kSx ~ 1 where 5x 
is the perturbation to the ion orbits in the fc-direction 
due to the electromagnetic fields of the instability. In 
another word, the magnetic field stops growing when the 
deviation of ion orbits is comparable to the instability 
wavelength. Here we give still another explanation of 
EqllOl as a condition that the current perturbation from 
the particle bouncing and trapping reaches the maximum 
possible value in one e-folding time. 

The Weibel instability is driven by the currents created 
by and located near the nulls of the perturbed magnetic 



field. This can be seen most simply by assuming that 
initially all particles move only in the ±y -direction with 
equal numbers moving at velocities —V y and +V y , re- 
sulting in zero j y everywhere. When there is a perturbed 
magnetic field B z with a k x , because of the sign change of 
B z near a certain null, the particles moving in one direc- 
tion are trapped (in the ir-direction) near the null and the 
particles moving in the opposite direction are deflected 
away from the null, creating a net current. Near the 
neighboring null, both the trapped and deflected parti- 
cles move in reversed directions, resulting a net current 
of an opposite sign. It is these alternating currents, rep- 
resented by the a-term in Eq[4l that drive the instability. 
The net current density j y created depends on how effec- 
tively the current of the opposite sign can be reduced by 
transverse deflection of particles away from the region of 
interest. The magnitude of the current can therefore be 
estimated from the deflected particle flux by 
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where jm — {\/2)n\qV y \ represents the maximum current 
density achievable, which would arise if all particles of 
one velocity sign leave the null region. 

Due to magnetic deflection, dv x /dt ~ uj c V y where uj c = 
qB z /Mc is the cyclotron frequency. Since dB z /dx « kB s 
near the null if B s is the magnetic field amplitude, we 
obtain 
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Comparing Eq. (jTUJ) and Eq. (fT2"]) , we can see that B s 
is the field strength such that the current growth in one 
e-folding time reaches ~ ju , after which the current can 
no longer increase and the mode saturates. This satu- 
ration mechanism also shows that q,V±, and M, in Eq. 
(flC)]) should be that of the locally anisotropic species (the 
ions in our case) since the net current creation from the 
trapping and bouncing would be zero for an isotropic 
distribution [29| . 

Combining Eq. (fTU|) with Eq. ([5]) or Eq. © we can derive 
the saturated level for the magnetic field B s for a small 
or large. In the small a limit, Eq. (fTT))) and Eq.© show 
that the largest B s occurs at a mode number of k m = 
fco / v5j which is smaller than that of the fastest-growing 
mode, fc*. The ratio of magnetic energy in this mode at 
saturation to the initial proton kinetic energy gives 
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Because of the mass ratio factor, in a single mode 
Es/Ekp ~ 2.6 x 10~ 16 for an electron-proton plasma 
with a — 1 and T e /m = Tix/M. This would also approx- 
imates the total magnetic field energy if the spectrum at 
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saturation is dominated by this mode. Because the mag- 
netic energy is only a small fraction of the proton kinetic 
energy and because the only way for electrons to gain 
energy is through the magnetic energy, only a negligible 
fraction of the proton energy can be expected to be trans- 
ferred to electrons when the Weibel instability saturates. 
The PIC simulations in Sec lIIII will indeed show that the 
electron energy gain is on the same order of magnitude 
as the magnetic field energy gain. The upshot is that for 
an initially unmagnetizcd ion-electron plasma with mod- 
erate anisotropy, the Weibel instability cannot generate 
electromagnetic fields of sufficient strength to energeti- 
cally couple them to the protons. 

For a ~ M/m, Eq.^ gives E B /E kp ~ 4.7 x 1(T 4 ~ 
0(m/M). To obtain more accurate results for large a, 
Eq.© and Eq. fTO]) can be used to obtain 

gs _ 1 (kc/uj pI ) 2 

E kp ~ 4 [1 + (kc/co pI ) 2 + M/m] 2 " 1 ' 

Eq. (jT4j) shows that in the large a limit the largest B s oc- 
curs at a mode number of = (u> p i /c) a/ M/m = uj pe /c, 
Es/Ekp ~ m/(16M). This result agrees with that 
from the calculation using the water-bag distribution [24|. 
Es/Ekp ~ m/(6M), within a numerical factor. 

In the literature, another saturation criterion was 
proposed[||, @, kp ~ 1, where p is the Larmor radius 
of one of the species in the plasma. For the Weibel 
instabilities driven by the electrons in a stationary ion 
back ground, this gives similar scaling for B s as EqfTOl if 
k = ujpe/c and p = p e , the electron Larmor radius, are 
used. However, for the proton-driven case here, Ref.[24| 
pointed out that this criterion with k = uj p i /c and p = pi, 
the ion Larmor radius, used in Ref.Q gave a larger B s 
than Eq .(PI). Recently, Ref. [6] claimed that this criterion 
gave the correct scaling, Es/E kp ~ m/M, if k = u> p i/c 
and p — p e were used. The analysis here shows that 
the largest B s occurs at k = u> pe /c, not k = u> p i/c. At 
k = Wpj/c, EqfT4l gives E B /E kp ~ (m/M) 2 . It is clear 
that Eq. (TlO|) is the correct saturation criterion for the 
magnetic field in the Weibel instabilities, valid for both 
electron-driven and ion-driven cases. 



III. PIC SIMULATIONS 

The saturation condition Eq. (|10[) was derived for a sin- 
gle mode, ignoring mode-mode interactions. To study 
how accurately Eq. (Ti"T)|) describes the saturation level 
when a range of unstable modes are present, whether 
Eq. (fl"3)) is a good estimation of the total magnetic field 
energy, and whether the electron gained energy is of the 
same order of magnitude as the magnetic field energy, we 
performed a series of PIC simulations with the fully ex- 
plicit PIC code OSIRIS [3(3]. The simulations also provide 
a check on the theory of Weibel instabilities in the low 
growth rate regime [261] and include any relativistic effects 
neglected in the analysis. This regime is motivated by 



turbulent astrophysical plasmas of interest around black 
holes which are expected to have low to moderate tem- 
perature anisotropy (a < M/m), moderately relativistic 
temperatures, and low growth rates. Previous PIC sim- 
ulations were in the large anisotropy and large growth 
rate (a > M/m and large £) regime@, |f3|. 

In practice, PIC simulations of low growth rate insta- 
bilities are more difficult than those of a high growth 
rate. The number of particles used in a PIC simula- 
tion plasma is usually much smaller than that of the 
actual plasma it simulates. The resulting inflated fluc- 
tuations in the simulations, including the associated col- 
lisions, set a practical lower limit on the growth rate of 
a collisionless instability that can be reliably measured. 
Eqs.([7]) and (TTSj) show that an actual proton-electron 
mass ratio of M/m = 1836 would produce too small a 
growth rate and too weak a magnetic field energy den- 
sity to be practically measured in a simulation. We have 
therefore used M/m = 10 — 30 in our simulations to 
check the ion mass scaling. Both electrons and ions were 
initially uniformly distributed in space, with a veloc- 
ity distribution of the form exp(— p x /p 2 x ) exp(— p 2 /p 2 v ) 
where p XtV is the normalized relativistic momentum com- 
ponents, p x , y = {v x , y /c)/^Jl - v 2 /c 2 - v 2 /c 2 . The z- 

component of the velocity was kept to zero. For the 
electrons, pt x — Pty — 0.41 were chosen, corresponding 
to ^/T e /mc 2 ~ 0.38 (T e = 74 keV) if we write T e /mc 2 w 

Ptx/\J l +Ptx + Pty For tne ions ' Ptx = 0.13 and p ty 
0. 1 1 were chosen, corresponding to */ Ti x /Mc 2 ~ 0.12, 
^jTjy/Mc 2 ~ 0.39, and a — 9. (For protons, this cor- 
responds to T Ix ~ 15MeV and T Jy ~ 145MeV.) The 
simulations used 1-1/2D (one full spatial dimension, x, 
and one half-dimension, y, where no spacial variation was 
allowed, and three velocity components but with V z re- 
maining zero throughout the runs) with a box size of 240 
c/Lu P e. The cell size was Ax = 0.1c/oj pe and time step 
was At = 0.099/wp e . To reduce fluctuations, 2500 to 
18000 particles per cell were used for each species in typ- 
ical runs. The total energy in these runs was conserved 
within 10~ 6 . 

Fig. Q] shows the Fourier spectrum of B z at two dif- 
ferent times, t\ = 166a->~ e 1 (in the linear growth stage) 
and ti = 305u p( } (after saturation), for the M/m = 10 
case (P = 3 and 2500 particles per cell were used.) Dur- 
ing the linear growth stage, the growth rate 7 for the 
k = 0Aui pe /c was measured to be 0.026w pe , agreeing rea- 
sonably well with the analytical result of 0.031w pe ob- 
tained directly from Eq. (Here, w 0.45 so the 
small ^/-expansion result of Eq. ([5]) is not accurate and 
the analytical result is obtained using exact values of the 
Z-function.) Following the same mode to saturation at 
t = 305/<x>j, e , the right hand side of Eq. p0[) was found to 
be 0.017ix>p e , about 65% of the linear growth rate, consis- 
tent with Eq. pop as a qualitative saturation criterion. 
This also implies that Eq. (flU)) is consistent only as an 
order-of-magnitude estimate. Fig. [T]also shows the shift 
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of the dominant mode to lower k as the instability satu- 
rated, since k* of the maximum growth mode decreases 
as a decreases and also k m of the highest saturation mode 
is smaller than k* . 

Fig. [2] shows the time evolution of the energy change 
in ions, electrons, and magnetic field (£? z -component), in 
units of the fraction of the initial ion energy. The curve 
for the magnetic energy is the total magnetic magnetic 
energy, with contributions integrated over the full range 
of k. This magnetic field energy reached ~ 1.8 x 10~ 3 
of the initial ion energy, which is of order 1.3 x 10~ 3 
estimated from Eq. (fT3"]) . The electron energy change 
at saturation was about 3.4 x 10~ 3 , significantly higher 
than the gain of 1.8 x 10~ 4 that corresponds to the gain 
from intrinsic collisions in the code: The latter baseline 
collisional gain was measured in a run with an isotropic 
ion distribution and found to depend linearly on t and 
inverse linearly on the number of particles per cell used. 
Therefore, the electron energy gain in the simulation of 
Fig. 2 was mainly due to collisionless processes as shown. 
After saturation, the electron energy remains at the same 
order of magnitude till the end of the simulation (t — 
3000^ or 948w~/). 

The ion mass scaling of (TIB"]) was checked by com- 
paring the M/m = 10 run with a run of M/m = 20. 
The saturated magnetic energy should scale as M~ 4 ac- 
cording to Eq. ([T3"l) , and the time to reach saturation 
should scale as M 3 / 2 (Eq. [7]). Because the ion-electron 
collisional energy exchange rate in our simulations was 
found to scale as M _1 , as in a real plasma [3l|, it was 
necessary to increase the number of particles per cell 
to 18000 in the M/m = 20 simulations to reduce the 
collisional energy exchange rate to below the collision- 
less instability-induced level. Fig. [3] shows the energy 
changes for the M/m = 20 case. The saturated mag- 
netic field energy reached at ~ 1.3 x 10~ 4 of the initial 
ion energy, somewhat less, but to within order of mag- 
nitude of the 5.3 x 10 -4 estimated from Eq. (fT3"|) (using 
(3 = 1.5). The electron energy change at saturation was 
~ 3.2 x 10~ 4 , roughly scaling as M~ 3 4 when compared 
with the M/m = 10 case. Overall, the PIC simulations 
support Eq. (| 1 3[) as an estimation of the saturated mag- 
netic field energy and that the electron energy gain is of 
the same order of magnitude as the magnetic field energy. 

IV. SUMMARY AND DISCUSSIONS 

In summary, the linear growth rate of the Weibel insta- 
bility driven by the proton temperature anisotropy has 
been calculated in the low anisotropy limit [Eqs.([5]) and 
([7])]. This growth rate is also used to derive an expression 
for the ratio of saturated magnetic field energy to the ini- 
tial proton kinetic energy, Eb /E^p, in the low anisotropy 
limit [Eq. fllj)) ]. These results show that Es/E^p scales 
as (m/M) 4 when ^jT e /m ~ yjT[ x /M and is extremely 
small for a — 1. 1.5D PIC simulations with a reduced ion 
to electron mass ratio support the scaling relations and 



also show that the electron energy gain from the proton 
anisotropy is comparable to the gain in magnetic field en- 
ergy. Although in three dimensions one additional mode 
polarization is allowed and the number of unstable modes 
can double, we do not expect orders of magnitude change 
in the electron energy gain. From these analyses and sim- 
ulations we conclude that the proton-driven Weibel insta- 
bility will not provide an efficient faster-than-Coulomb 
electron-ion temperature equilibration mechanism in an 
intially unmagnetized two temperature plasma with low 
to moderate anisotropy. 

We have also shown that the low anisotropy formulae 
derived here approach to the previous high anisotropy 
results when a — > M/m [Eqs.© and l|14p]. This makes 
explicit the implied assumption of a > M/m used in 
the discussion of magnetic field generation in collision- 
less shocks through the Weibel instability [f| |24|. It 
also provides guidance on how to properly scale PIC 
simulations results with reduced M/m. For example, 
PIC simulations on shocks in galaxy clusters with Mach 
number M s h — 20, electron thermal velocity Vth,e — 
0.05c, and ion thermal velocity Vth,e = 0.005c were 
performed using M/m = 100 in Ref.[6j. Since a — 
Ms h (m/M)(v t h,e/v t h,i) 2 = 400 > M/m, the simulations 
showed an Es/E^p ~ 10 -3 , consistent with m/(16M) 
predicted from Eq. pT)) . However, if the same M,v t h,e> 
and vth,i were used for M/m = 1836, a — 22 < M/m 
and Eq.([T3"|) would predict an E B /E kp ps 4 x 10 -6 . 

Finally we wish to point out that the situation may 
be greatly different for an initially magnetized plasma 
in which electrons are confined to local magnetic field 
lines transversely. For modes with k perpendicular to 
the equilibrium Bo, the electrons would be practically 
immobile in the direction perpendicular to both Bo and 
k. As explained in SecJITl the Weibel instability is 
driven by the current induced by the perturbed mag- 
netic field B z , represented by the term proportional to 
q in Eqf?] However, additional currents with a differ- 
ent phase are also driven by the accompanied electric 
field E y , which are represented by the second and third 
terms on the left of Eqf?] These currents are present even 
for an isotropic distribution and oppose mode growth. 
If the electrons were magnetized and could not respond 
to E y to produce a curent, the second term would not 
be present and the growth rate would be independent 
of m/M. This indicates that the required equilibrium 
magnetic field should be such that uj ce >> uj p j >> u> c j 
for this to happen, where to ce and u> c i are the electron 
and ion cyclotron frequencies in the equilibrium mag- 
netic field. In the case of the ion-driven (electricstatic) 
two-stream instability with magnetized electrons, the re- 
sults were indeed found to be independent of the mass 
ratio [12 • If this also were true for the Weibel instability 
more generally, the saturated magnetic fields could be 
E B /E kp re 8 x 10- 3 cv 5 /[(i + a) 5 (3 + 2a)] = 2.5 x 10" 4 
for a = 5.4, significantly larger than that predicted by 
Eq. lfTB")) . However, the electrons still cannot gain signifi- 
cant energy even in this enhanced magnetic field. 
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FIG. 2: The energy change of the ions, electrons, and mag- 
netic field, in unit of the fraction of the initial ion energy, for 
the M/m = 10 case. 
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FIG. 3: The energy change of the ions, electrons, and mag- 
netic field, in unit of the fraction of the initial ion energy, for 
the M/m = 20 case. 



